Field-dependent thermal transport in the Haldane chain compound NENP 
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We present a study of the magnetic field-dependent thermal transport in the spin S = 1 chain 
material Ni(C2H8N2)2N02(C104) (NENP). The measured thermal conductivity is found to be very 
sensitive to the field-induced changes in the spin excitation spectrum. The magnetic contribution to 
the total heat conductivity is analyzed in terms of a quasiparticle model, and we obtain a tempera- 
ture and momentum independent mean free path. This implies that the motion of quasiparticles is 
effectively three dimensional despite the tiny interchain coupling. 

PACS numbers: 75.40.Gb 66.70.+f, 75.47.-m 



The recent theoretical interest concerning heat trans- 
port in one dimensional (ID) spin systems (see review ar- 
ticle [1] and the most recent papers [2, 3, 4, 5, 6, 7, 8]) was 
greatly stimulated by observations of a large spin ther- 
mal conductivity k s in two-leg Heisenberg 5 = 1/2 ladder 
compounds (La,Sr,Ca)i4Cu2404i at high temperatures, 
with the mean free path of spin excitations I reaching 
3000 A [9, 10, 11]. However, experiments on several 
Heisenberg 5 = 1 chain compounds [12, 13, 14] pro- 
vided evidence for considerably lower n s , with I < 60 A. 
This is surprising because the 5 = 1 chain model and 
the 5 = 1/2 ladder model are essentially equivalent [15]. 
Both adopt a spin liquid state with exponentially decay- 
ing correlations and an energy gap in the spin excitation 
spectrum. 

The spin thermal conductivity k s at a finite frequency 
lo can be decomposed into a singular and a regular part 



ReK s (uj) = D th S(uj) + K reg (w), 



(1) 



where D t h is the thermal Drude weight. In our exper- 
iment, we measure the dc conductivity Re /« s (w = 0). 
For integrable models and continuum field theories, heat 
transport is ballistic even at finite temperature T, D t h > 
0, as conservation laws prohibit the decay of the heat 
current. Both extrinsic sources of scattering, such as de- 
fects and phonons, and intrinsic spin-spin interactions 
which are not integrable, render the heat conductivity 
finite [1, 2, 3, 4, 5, 6, 7, 16, 17]. 

Measurements of the thermal conductivity in external 
magnetic fields, which can strongly modify the spin exci- 
tation spectrum, offer detailed information on scattering 
mechanisms limiting k s . However, magnetic fields typi- 
cally available are too weak to noticeably influence the 
spectrum of spin excitations in the previously investi- 
gated [12, 13, 14] S = 1 chain materials AgVP2Sg and 
Y2BaNi05 with strong intrachain exchange J > 250 K. 

In this Letter, we present results on field-dependent 
heat transport in one of the model low- J 5 — 1 chain ma- 
terials, Ni(C 2 H 8 N 2 )2N02(G104), viz. NENP. The mean 
free path of the spin excitations, evaluated from our data, 



is large and temperature- independent allowing us to iden- 
tify the most relevant scattering mechanism. The heat 
transport at low temperatures is determined by rare de- 
fects, cutting the spin chains into segments, and not by 
the intrinsic interactions. 

NENP crystallizes in the orthorhombic Pnma space 
group with lattice parameters a=15.223 A, 6=10.300 A, 
and c=8.295 A [18]. The 5 = 1 spins of Ni 2 + form 
chains along the 6 axis with exchange constant J « 43 K 
[19], while the interaction ,/' between the chains is much 
weaker (J'/ J ~8x 10~ 4 according to [20]). Therefore, 
low-temperature 3D ordering is neither expected nor ob- 
served. Neglecting the interchain interaction for the mo- 
ment, the appropriate Hamiltonian is 

H = ]T{JS 4 S m + D{Sl) 2 + E[{Slf - (S l y ) 2 } 

i 

+ A^-gB}, (2) 

where D and E are single- ion anisotropy constants, B 
is the magnetic field, and g is the gyromagnetic tensor. 
In an ideal isotropic antiferromagnetic (AFM) 5 = 1 
chain, the excitations are triply degenerate with a gap 
A = 0.41J at the AFM wavevector &af = n/d, where 
d = 6/2 is the distance between neighboring spins along 
the chain. In NENP, the strong planar anisotropy and 
weak orthorhombic anisotropy (D/J = 0.2 and E/Ja 
0.01 [18, 21]) split A into three gaps A, (i = 1,2,3) 
with the zero-field values A? ps 29 K, A^ = 14.3 K, and 
A§ = 12.2 K [20]. With increasing B || 6, Ai stays con- 
stant, A 2 increases, and A3 decreases such that it should 
close at the critical field B c « 10 T and the system should 
enter a gapless Luttinger liquid (LL) state for B > B c 
[22]. This does not happen, however, because the chemi- 
cal environment of every second Ni atoms along the chain 
is oriented in a different direction [23]. This alternating 
tilting introduces an additional term J2i(~ -g s tB 
in the Hamiltonian Eq. (2), where the staggered trans- 
verse field B st = I gstB I is proportional to the homoge- 
neous field B. As a consequence, the gap remains finite 
at B = B r and increases above B r . 
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FIG. 1: (a) Thermal conductivity of NENP parallel and 
perpendicular to the spin chains as a function of temperature 
at B = and B = 8.2 T. (b) The relative change of the 
thermal conductivity of NENP parallel to the spin chains as 
a function of magnetic field at several constant temperatures. 



The crystals of NENP used for our experiments were 
grown as described in Ref. [24]. One crystal of dimen- 
sions 1.1x3x0.5 mm 3 with the longest dimension along 
the b axis was used for the measurements of the thermal 
conductivity along the spin chains. A sample of dimen- 
sions 1.8x0.8x1.4 mm 3 for measurements perpendicular 
to the chains was cut from another crystal. For the ther- 
mal conductivity measurements, we employed a standard 
steady state method with the same arrangement of the 
thermometers and heater as described in Ref. [25]. The 
experiments were performed in the temperature range be- 
tween 0.3 and 50 K in magnetic fields up to 16 T applied 
parallel to the chain direction. 

The thermal conductivity k(T) both parallel (k") and 
perpendicular (k^) to the chain direction in zero field and 
in B = 8.2 T are shown in Fig. 1(a). The relative changes 
of k" as a function of magnetic field k^(B)/k^(0) at sev- 
eral constant temperatures are displayed in Fig. 1(b). 
The striking observation is that a magnetic field leads 
to a strong enhancement of the thermal conductivity up 
to 5 times its zero-field value. This strong enhancement 
is restricted to the direction parallel to the chains. The 
small increase of k _l (-B) can easily be attributed to a less 
than 1° deviation of the heat flow direction from being 
exactly perpendicular to the b axis. 

The total thermal conductivity of a magnetic insula- 
tor can be represented as k = K p h + k Si where the two 
terms on the right-hand side correspond, respectively, 
to the phononic and magnetic contributions to the heat 
transport (possible spin-phonon drag contributions are 
included in k s ). The spin excitation spectrum at B = 
is gapped with the smallest gap A3 = 12.2 K, there- 
fore at T <C A3, both spin thermal conductivity and 
phonon-spin scattering are negligible. With increasing 
B at a constant temperature, A 3 (B) decreases and the 
number of thermally activated spin excitations increases. 
Therefore, k s (B) should increase, while Kph 

(B) should 

decrease because of the growing phonon scattering by 
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FIG. 2: Magnetic contribution to the thermal conductivity 
of NENP parallel to the spin chains as a function of temper- 
ature at several constant fields. The solid lines represent the 
calculated k s (see text). The error bars arise from the uncer- 
tainty in subtraction of K p h from the total measured thermal 
conductivity. 



spin excitations. Because of the quasi-lD nature of the 
spin system in NENP, njr is small and «; « K^ h . Thus, 
the observed negligible influence of the magnetic field on 
k 1 - suggests that spin-phonon scattering is weak. The 
increase of k^(B) in fields < B < B c clearly demon- 
strates that all field-induced changes in k" originate from 

A salient feature of the curves shown in Fig. 1(b) is 
that for all temperatures below about 1.5 K there is a 
region of low fields where k"(B)//«"(0) are practically 
field-independent. This means that below about 1.5 K 
k" (B = 0) is purely phononic. Therefore, by subtracting 
the zero-field values from k" at these temperatures we 
obtain a good estimate of k s (B). 

We have measured k"(T) below 1.5 K at several con- 
stant fields. The spin contribution k s (B, T) = k" (B, T) — 
k" (0, T) is shown in Fig. 2. For a quantitative analysis of 
the data, we consider the heat transport associated with 
excitations from the singlet ground state to the lowest 
triplet branch. As the energy gaps for the other two 
branches of the triplet either increase with B or stay 
constant, their contribution to the heat transport below 
T w 1.5 K can be disregarded. The interaction between 
chains is very weak in comparison with the intrachain in- 
teraction; nevertheless, it still leads to a dispersion per- 
pendicular to the b axis with a bandwidth of about 2 K 
and 0.8 K along the a and c axes, respectively [20, 26]. 
The dispersion we use in order to analyze k s in low fields 
(hence the subscript "If") lT<_B<6Tis given by 

£l/ (k) = {[(A°+A°)/2] 2 + V 2 (M-7r) 2 
+ (A£ Q ) 2 [(l + cos£; a a)/2] 2 
+ (A£ c ) 2 [(l + cosfc c c)/2] 2 }^ 2 - 56MB B, (3) 

with constants A^ = 14.3 K, A° 3 = 12.2 K, AE a = 7.5 K, 
and AE C — 5.0 K taken from neutron scattering exper- 
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FIG. 3: Mean free path of spin excitations, calculated for 
1 T < B < 6 T. The error bars correspond to those of k 3 
shown in Fig. 2. An approximately constant mean free path 
can be obtained from a model (solid line) where strong impu- 
rities effectively cut the spin chains into segments of length 
I. In contrast, the data cannot be described (dashed line) by 
weak fluctuations of J. 



iments [20, 26] and g b = 2.15 from an ESR study [19]. 
The value of V = 2.49J, used in the analysis, is pre- 
dicted by theory [27] and is confirmed by the inelastic 
neutron scattering measurements [20, 28]. Equation (3) 
describes correctly the linear-in-i? decrease of the energy 
gap at q — (0, Tr/d, 0) between 1 and 6 T observed in 
the ESR experiment [19], but does not account for de- 
viations at higher fields caused by the influence of the 
aforementioned transverse staggered field. 

Within the Boltzmann equation approximation, the 
spin thermal conductivity parallel to the chain direction 
is represented by 

7r/a 7r/c Tt/b 


where n = 4 is the number of spins in the unit cell of 
NENP, / is the distribution function, e(k) is the energy, 
and I is the mean free path of spin excitations. In Eq. (4), 
we use the Fermi distribution / = (exp(e/fcsT) + 1) _1 
(taking into account the hard-core repulsion of the ID 
excitations) [29], but as the temperature is well below 
the gap, one can equally well use a Boltzmann- or Bose 
distribution. 

In the field region 1 T < B < 6 T, where e(k) = ej/(k), 
the only unknown parameter in Eq. (4) is I. Thus, the 
mean free path, 1{B,T), as a function of B and T can 
be extracted from the experimental data of Fig. 2 and is 
shown in Fig. 3. Within the experimental accuracy, the 
mean free path is independent of both T and B, with the 
average value I = 0.75±0.1 /im. Remarkably, / for NENP 
is as large as the highest values of the mean free path 
found for S = 1/2 Heisenberg chains and ladders [10, 25, 
30], where impurities are the main source of scattering 
for spin excitations at low T. Both for spin-spin and 



spin-phonon scattering, the mean free path is expected 
to increase rapidly with decreasing T [3]. We therefore 
conclude that, at least at T < 0.04 J, scattering by defects 
and not the intrinsic interactions determine k s . A mean 
free path due to the intrinsic processes is theoretically 
expected to be huge [2] as Umklapp processes which relax 
momentum are exponentially supressed for T< J. 

A field and temperature independent mean free path 
cannot be obtained from a purely ID model. For both 
weak and strong impurities, one expects [25] in one di- 
mension a momentum dependent mean free path propor- 
tional to the square of the velocity, Ik oc uf, implying 
effectively a variation of / linear with T. The situation 
is, however, completely different when one takes the tiny 
three dimensional coupling between the chains into ac- 
count. If spin excitations can pass "strong" impurities, 
which effectively cut the chains into segments, by hop- 
ping to the next chain instead of tunneling through the 
defect, one naturally obtains a mean free path given by 
the distance of the defects. We have calculated the weak 
T dependence of I (solid line in Fig. 3) assuming a small 
density of local, infinitely strong potential scatterers in 
a model defined by the dispersion relation (3) for a gap 
A of 5 K (for T <C A the result is almost independent of 
A). Such a simple calculation to linear order in the den- 
sity of defects is valid as the interchain coupling J 1 (or 
more precisly the bandwidth in perpendicular direction) 
is much larger than the inverse of the time needed by a 
spin excitation with energy T to propagate to the next 
defect, J' > VTJd/l. 

The calculations for the mean free path dominated by 
strong impurities describe the data very well (solid line in 
Fig. 3). This is not the case when one assumes that small 
fluctuations 5J < T of J dominate transport (dashed 
line in Fig. 3). In the latter case, the scattering rate is 
proportional to the density of state [obtained again from 
Eq. (3)]. 

The model we used for the low-field calculations fails 
to describe k s (B,T) at B > 6 T. This is illustrated in 
Fig. 4(a), where the dashed line is n s (B) calculated for 
a constant temperature T <C J using Eqs. (3) and (4) 
with 1 = 1. The calculated k s (B) shows an increase when 
B approaches the LL state from the gapped state, with 
a plateaulike feature at B c broadened by the interchain 
interaction. This is the expected generic behavior for a 
ID spin system near a quantum phase transition from a 
gapped to a gapless state and has, indeed, been observed 
for the S = 1/2 chain compound CuPzN in [25]. The dif- 
ferent behavior of NENP for high fields arises obviously 
from the energy gap induced by the staggered field. 

We have fitted Eq. (4) to the k s (T) data presented 
in Fig. 2 using a modified dispersion relation e(k) — 

iJAj_ + (k) , where £//(k) is given by Eq. (3) and 

the energy gap Aj_ is induced by the staggered field. 
A similar heuristic fitting formula has e.g. been used 



4 



3 | = — | 20 




B (T) B (T) 

FIG. 4: (a) K a (B) calculated for T = 1.33 K with (solid 
line) and without (dashed line) consideration of the staggered 
transverse field. The open circles represent the experimental 
data, (b) The open squares are the minimum energy gap 
A m i n estimated in our analysis. For comparison, the energy 
gap obtained from specific heat measurements from Ref. [32] 
are show (stars). The solid squares are the fitted values for 
Ax (the vanishing of Ax below 6T may be an artifact of the 
fitting procedure as A m i n depends very litte on Ax in this 
regime). Inset: The fit values of the mean free path of spin 
excitations. 



in Ref. [31] to describe the S = 1/2 chain compound 
CuCl 2 -2((CD 3 ) 2 SO) where, similar to NENP, a staggered 
field induced by a uniform magnetic field leads to a finite 
gap in the spin excitation spectrum. For each value of B 
between 1 T and 16 T, Eq. (4) was fitted to the exper- 
imental k s (T) data, shown in Fig. 2, with two free pa- 
rameters I and Ax , assuming that I is T- independent for 
each B. The resulting k s (T) curves are shown in Fig. 2. 
The fit values for A±(B) are shown by solid squares in 
Fig. 4(b), the data for 1(B) are presented in the inset of 
Fig. 4(b). The solid line in Fig. 4(a) is calculated using 
these data for 1(B) and A±(B). As shown in Fig. 4(b), 
there is an agreement between our data for the minimum 
energy gap A min = [((A« + A°)/2 - g b ^ B Bf + A 2 ± ]V2 
at k = (ir/a,7r/d,n/c) and the data for the energy gap 
obtained from the specific heat measurements [32]. The 
essential result is that the mean free path I remains close 
to its low-field value I — 0.75 fim in the entire 1-16 T 
field region. This is again consistent with the notion of 
a mean free path determined by strong impurities. 

In summary, from our measurements of the anisotropic 
thermal conductivity of the S = 1 Haldane chain com- 
pound NENP, we identify a large magnetic contribution 
along the spin chain direction. The mean free path of 
spin excitations is orders of magnitude larger than previ- 
ously observed for other S — 1 chain materials and of the 
same order of magnitude as in the best 5=1/2 chain 
and ladder compounds. We have argued that the absence 
of a temperature and field dependence of the mean free 
path can be explained by rare defects, which effectively 
cut the spin chains into segments, in combination with 
a tiny interchain coupling. The measured values of spin 



thermal conductivity may also serve as a lower limit for 
future theoretical estimates of the intrinsic diffusive con- 
tribution to the heat transport in S = 1 AFM chains at 
low temperatures. 
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